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ABSTRACT 

We  consider  the  steady  two-dimensional  flow  under  gravity  of  water  from 
one  reservoir  to  another  reservoir  through  a porous  rectangular  isotropic 
homogeneous  dam  with  impervious  bottom.  Using  the  maximum  principle  we 
give  a proof  of  the  convexity  of  the  free  boundary. 
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SIGNIFICANCE  AND  EXPLANATION 


We  consider  the  steady  two-dimensional  flow  under  gravity  of  water  from 


one  reservoir  (on  the  left)  to  a lower  reservoir  (on  the  right)  through  a 


porous  rectangular  isotropic  homogeneous  dam  with  impervious  bottom 


Because  of  gravity  the  water  does  not  flow  through  the  entire  dam  and  the  dam 


is  dry  near  its  upper  right  corner.  The  interface  separating  the  dry  and 


wet  regions  of  the  dam  is  a free  boundary.  Recently,  Friedman  and  Jenkins 


have  proved  that  the  free  boundary  is  convex.  We  give  a different  proof 


which  uses  only  the  maximum  principle  and  its  generalizations 
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A PROOF  OF  THE  CONVEXITY  OF  THE  FREE  BOUNDARY  FOR  POROUS  FLOW  THROUGH 


A RECTANGUIJVR  DAM  USING  THE  MAXIMUM  PRINCIPLE 

* 

C.  W.  Cryer 


1 . The  p«un  Problem 

We  consider  the  following  Dam  Problem.  Water  flows  steadily  under  gravity  through 
a porous  rectangular  isotropic  homogeneous  dam  ABCF  with  impervious  bottom  BC  of 
length  L from  a reservoir  of  height  H to  a reservoir  of  height  h (see  Figure  1.1). 
The  water-air  interface  AE  is  a free  boundary  which  we  denote  by  T. 


Figure  1.1  : PI ow  through  a r ectangular  dam . 

The  mathematical  problem  is  as  follows  (Bear  (1972),  Baiocchi  (1972), 
Baiocchi  and  Capelo  (1978))  ■ 

Find  functions  ^(x)  (the  height  of  the  free  boundary)  and  u(x,y)  (the 
hydraulic  head)  such  that  (from  the  equation  of  continuity  and  Darcy's  law) s 
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I—  - u * u = V'u  = 0,  in  fl  , 

pyj  xx  yy 

(1.1) 

together  with  the  boundary 

condition • 

• 

u - H, 

on 

AB, 

(interface  with  water  at  rest)  , 

(1.2) 

S"0* 

on 

BC, 

(impervious  boundary)  , 

(1.3) 

u « h. 

on 

CD, 

(interface  with  water  at  rest)  , 

(1.4) 

u - y, 

on 

DE, 

(interface  with  air)  , 

(1.5) 

u » y. 

on 

EA, 

(interface  with  air)  , 

(1.6) 

on 

EA. 

(streamline)  . 

(1.7) 

Here,  (1  is  the  (unknown)  domain, 

S)  ■ { (x,y)  : 0 < y < <fi  (x)  , 0 < x < Ll  , 
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and  denotes  the  unit  outward  normal  derivative.  The  physical  significance  of  u 


is  that 


u - y + p/pg  , (1.9) 

where  q is  the  acceleration  due  to  gravity,  p is  the  fluid  pressure,  and  p is  the 
fluid  density. 

The  dam  problem  is  a well-known  problem  in  the  theory  of  porous  flow.  Cryer  11976, 
p.  54)  summarizes  the  numerous  numerical  solutions.  There  are  also  three  analytical 
solutions  involving  elliptic  integrals:  the  first,  is  due  to  Davison  (1932,  1936, 

1936a|:  the  second  is  due  to  Hamel  11934)  and  is  described  by  Muskat  11937,  p.  3031 
and  Bear  (1972,  p.  398):  the  third  is  due  to  Polubarinova-Kochina  (1962,  p.  284) 
(concerning  misprints  see  Cryer  (1976,  p.  54)).  Although  these  analytical  solutions 
are  rather  complicated,  they  can,  and  have,  been  evaluated  numerically. 

Despite  the  fact  that  the  analytical  solution  of  the  dam  problem  is  known,  the 
problem  is  still  frequently  considered  in  the  literature  because  it  serves  as  a useful 
model  problem  for  porous  flow  free  boundary  problems.  In  a pioneering  paper, 

Baiocchi  (19721  reformulated  the  problem  as  a variational  inequality  for  the  function 


and  derived  many  properties  of  w.  Further  properties  of  w and  u • y - w have 


been  obtained  by  Caff are  111  I to  appear),  Friedman  I 197bl , Jensen  11977),  and  Friedman 


and  Jensen  11977).  Aitchison  11972]  gives  an  expansion  for  the  solution  near  the  separa 


tion  point  E,  and  Aitchison  11977)  gives  numerical  solutions 


Recently,  Friedman  and  Jensen  11977)  have  proved  that  f is  convex  under  the 


following  assumptions  (which  had  previously  been  proved  in  the  above-mentioned  work  of 


Friedman,  Jensen,  and  Caffarelli) 


is  strictly  monotone  decreasing  for  0 < x < L 


u is  bounded  on  il 


Making  the  same  assumptions  we  give  a different  proof  of  convexity  which  uses  only  the 


maximum  principle  and  its  generalizations.  Some  preliminary  results  of  general  appli 


viability  are  given  in  Section  2,  and  are  then  applied  to  the  case  in  hand  in  Section  3 


The  basic  ideas  in  Section  3 have  been  known  to  us  for  almost  twenty  years  and  were 


originally  motivated  by  proofs  of  convexity  for  fluid  dynamics  free  boundary  problems 


(Birkhoff  and  Zarantonello  11957,  p.  84),  Gilbarg  (1980,  p.  373)1.  Since  the  approach 


is  effective  in  many  fluid  mechanics  free  boundary  problems  we  hope  that  the  same  will 


be  true  for  porous  flow  free  boundary  problems 


The  present  paper  illustrates  the  power  of  the  maximum  principle  as  a tool  for 


analyzing  free  boundary  problems.  It  is  appropriate  to  men' ion  that  the  maximum 


principle  was  probably  first,  applied  to  the  dam  problem  by  Davison  11938)  and  Shaw 


and  Southwi  II  (19411.  Some  additional  references  to  the  application  of  the  maximum 


principle  to  free  boundary  problems  are  given  by  Cryer  (1977,  Section  11.13) 


I 


. I'leliminaty  Kesu  1 1 s on  the  Maximum  It  m.  ij>lf  and  I'  i f f or out  i.«  1 Geometry 

We  will  requite  nrveul  vt'isinns  of  l In-  maximum  | i in  lple  foi  a 11’  a 1 function 
u « Cdl)  which  is  harmonic  in  a bounded  domain  1.  with  boundary  30  in  the  complex 
z * x * ty  - plane.  Wc  assume  that  u is  continuous  m 0 except  perhaps  at  a finite 

number  of  loints  z ,z  t ,10. 

1 m 

The  Hoyt  Principle  il'rottet  and  Wei nbeiger  Jl‘tt-7,  p.  t»*>|  , Gilharg  and  Trudinqet  11977, 
p.  ill)  : l.et  u(z)  M in  0,  and  ulz^l  “ M at  a point  z ( r 30,  the  boundary  of 

0.  let  lil  u be  continuous  at  z and  Hi)  there  exist  a ball  !»  C i;  with  z t IB. 

0 0 

Then  the  outer  normal  derivative  of  u at  z . if  it  exists,  satisfies 

~ (z  ) ' 0 . Cl 

3n  0 

The  Maximum  Principle  states  that  if  u « C(0)  and  u y M on  10  then  u y M 
in  0.  In  the  present  context,  this  principle  is  not  adequate  because  we  sometimes 
consider  harmonic  functions  which  are  jxissibly  not  cont  inuous  at  the  lx  Hindary  ixMiits 
* A,  z,  * P,  and  z(  • K (Figure  1.1).  We  will  say  that  u satisfies  the 
Generali zed  Maximum  Principle  if 


impl les 


Kemai  k . 1 


u(«)  , N,  z*30j  z ^ z ^ , 1 < j < » , 


11  ( Z 1 x M 1 ll  0 . 


If  vIO  is  harmonic  in  the  unit  circle  p and  v is  continuous  on  P except 
at  C • 1,  it  is  no?  necessarily  true  that  v satisfies  the  Generalized  Maximum 
Principle  as  is  shown  by  the  following  examples , 


v(t)  ■ v(ro  ) • v(r,t')  * Real 


L-Li]  . » - r2) 

1 ' ll  I * rJ  - Jr  cos 


»•»  .i  function  which  is  dft'iniHl  on  the  unit  circle  P ami  v ■ 0 on  dD  except  at 


the  point  f.  ■ ♦!,  but  v(r,0>  * as  r * l. 
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V ( Cl  * v ( t . <*  > 


1 mag 


Jr  U - r Uln  8 


2 2 ' 

1 1 ♦ r - 2i  con  0 1 

satisfies  vli.tM  *0  as  r •*  1 for  alt  <) . O 

In  order  (or  the  Generalised  Maximum  Principle  to  apply,  u must  satisfy  sonw 
boundedness  condition.  In  articular,  it  is  known  (Goluzin  tided,  p.  2t>7|)  that  it  u 
is  bounded  above  in  il  then  u satisfies  the  Generalized  Maximum  Principle.  The 
follow  mu  lemma  uives  another  condition, 
l.ommu  2.  I 

let  v be  the  harmonic  conjugate  of  u in  ft.  (v  is  well-defined  and  single- 

valued  since,  by  assumption,  i!  is  a domain  and  hence  simply  connected.!  Let  3i!  be 

a rectifiable  closed  Jordan  curve.  If  I v I < M in  il  for  some  constant  M then 

1 v v 

u satisfies  the  Generalized  Maximum  Principle. 

Proof:  bet  * * «i(C)  denote  the  conformal  mapping  of  the  unit  disk  P * ( C t U|  < It 

m the  C-plane  onto  si.  The  functions  u(z)  and  v(zl  are  harmonic  conjugates  m si, 
and  so  the  functions  u(m(()l  and  viu>(Cl)  are  harmonic  conjugates  in  P. 

We  tecall  (Goluzin  11969,  p.  JHkll  that  the  Hardy* t ebesgue  space  h,  consists  of 
the  class  of  functions  f defined  on  P,  which  are  harmonic  on  P,  and  which  satisfy 

2« 

/ | ♦ (r . 6) | *d8  < m 

0 


for  some  constant  m and  all  re  (0,1). 

Since  v is  bounded,  v ho ( C.  1 1 e h.j  hence,  by  the  Riesz  theorem,  uW  ld)  « h, 
(Goluzin  ( 1969 , p.  192)1. 

bet  ♦(r.O)  * 4(0  ■ u(w(d).  Since  4 f h,  the  boundary  values  ♦(1,8)  ate 
defined  almost  everywhere  and  Poisson’s  formula  holds  (Coluzin  not'd,  p.  '>'1,  Theorem  ' 1 . 
Kudin  1 1 •>«.». , p.  232)1;  that  is, 


bet  u satisfy  iJ.l).  Since 
associates  sets  of  measure  zero  on 
p.  420,  Theorem  . Thus  | ♦ ( 1 . 0 ) I 


2"  2 

/ ♦11.0) — dt»  . 

0 \ - Jr  contit  ♦ r 

,V  \ «?  ,i  i tv* t\t  \ ab l e Jonhin  . ntvo,  the 

with  sets  of  nu'.isiu  e eero  on  }P  looluctn 

v M a.e.  on  .'P.  Kenvmhet  uw  that  the  weight 

-s- 


1 viUi  t ion  i •'■•****''  * * foimula  ih  non*negal I ve t it  follow*  fi\>m  I'o i »*v>n * * formula 


alH'vr  that  o M on  P,  i 

Koma  i V 

If  u anil  v a»e  haimonic  coni uoatox  and  v i*  hounded  then  u need  not  be 
bounded  a*  whown  by  t ho  example 


v»  ♦ i v - i t n ( x ♦ \ y ) . 0 

Furthoi  'nmrnt*  on  t ho  enoialisod  Maximum  Principle  will  he  found  in  the  appendix. 

In  the  dan  pioblem  w«*  ai  e given  the  value*  ot  u and  u on  Tho  following 

tv 

lemma  tunun.il  irrs  the  i o 1 at vonsh l p?  between  tho  der i vat ive*  of  a funct  ion  on  .t  cuive, 
iht*  lemma  diffoi*  fr»vn  tho  u.tti.tl  Frenet  formula.**  (Ftnenhart  U'^v\  p,  „'•*)>  in  that 
tho  . utv.itmo  i may  bo  eithei  positive  oi  negative. 

Lemma  - . 2 


lot  tho  boundary  It'  of  i?  bo  defined  par amet  r ioally  by  x - ti*',  y - n(*', 
whoio  * denote*  ate  length  along  in  tho  lenitive  dnection  *o  that  >}  i * to  tho 

left  of  .V  (Figure  .'.l'. 

lot  ^ • d(*  d*  • v denote  the  xigned  ou feature,  *o  that  Jf  *.  »*  than  >J 

i*  convex.  Let  t be  the  unit  tangent,  and  n the  unit  outward  normal,  riven 

t - (t  ^ ,t  - (♦*.  , MV>  , 

n • (n^  .n  J - (♦n,-(,)  ♦ 

t*  ♦ h*  • 1 • it  ♦ »\n  - o , 

* • Cn  - tn  , 

i • -*n*  u • ♦ tri  , 


d* 


dn 

ds 


- *vt  , 


lot  c - »f(x,y>  be  twice  oont  l nuoual  y differentiable*  We  denote  by  the 


derivative  of  %*  along  ' : that  is. 


vV  (£  i s > , n (a)  > 

dn 


Wo  denote  bv  c and  W*  tho  duoctional  de»  ivat  tve*  ot 
n t 


along  n and  t*  that  is. 


ssi  i i^irrr 


1,  Convexity  .’t  the  Free  Boundary 

The  basic  idea  is  to  cons  1 de  l the  function 

U - U - V i . 


u.  and  then  derivatives,  all  of  which  are  harmonu  functions 

The  solution  u can  be  reflected  by  symmetry  across  IV  (Courant  and  Hilbert 
IWt:.  V'.  3721)  so  that  u is  infinitely  differentiable  neat  PC.  It  thus  follow, 
from  the  basic  regularity  theory  for  elliptic  equations  (O.ilbaig  and  Trudmoei  jl-v'. 
Section  e . 4 ] 1 that  u vs  twice  continuously  differentiable  on  except  possibly  at 

F . In  this  section,  statements  alwut  the  values 
of  the  derivatives  of  u and  U on  3il  should  be  undet stood  to  exclude  tie  joints 
A,  D,  and  E. 

Using  (1.3)  through  (1.7)  direct  comput at  ion  yields: 


the  points  2^  - A.  2,  « P,  and  2? 


,'n  FA : u * ♦V'(x)  [1  ♦ ic'ix))") 


•W’lx))*  ’ll  * (*• ' (x) ) ‘ t 
I gi  ad  III*  - 1 g . 


On  ABCl):  u * 0,  C*  * grad  1’!'  ' 1 t . 


■■'ll  DE : u - 1 , 
V 


* | grad  V " '14. 

The  fact  that  0 is  constant  on  1'  means  that  m porous  flow  v'  plays  a 
which  is  similar  to  that  played  by  the  velocity  g in  fluid  dynamics. 


l'.4' 


Since  it  was  assumed  1(1.14''  that  is  bounded,  we  ma\  appl)1  the  Oenet  al  i red 


Maximum  Principle  to  u^  and  obtain 


y 


1 . in 


Now  consider  u.  Since  u < C(il'  the  Maximum  Principle  implies  that  u attains 
its  extrema  on  .hi,  but,  by  the  Hopf  Principle,  these  extrema  cannot  be  attained  on 
BO  because  ,lu  .'n  * 0 there.  It  has  been  assumed  that  is  nenct  oni.  al  ly  decreasing. 


sc>  that  h s u ' II  on  P.  It  follows  that  h ' u ' II  on  CPFAB  and  so 

h ' u ' H in  •'  . 


i ) .e' 


i 

I 


I 


». 


y- 


Next,  we  consider  ux-  By  the  Hopf  Principle  applied  to  u on  AB  and  CD  we 
See  that 


u • 0 on  AB  U CD 


(3.7) 


Where  we  have  used  the  fact  that  u can  be  reflected  across  BC.  The  boundary  conditions 
for  u together  with  (3.b)  imply  that  the  harmonic  function  u - y is  non-negative 
on  But  u " V * 0 on  DE  and  hence,  by  the  Hopf  Principle, 

ux  v 0,  on  DP  . (3.8) 

On  BO  we  have  u 
That  is. 


e u = 0;  noting  (3.5),  the  Hopf  principle  implies  that.  (u  ) N 0 

> y y - 


” (u  ) < 0,  on  BC  . 

* * YY 

But  u is  smooth  on  BO,  u^iB)  « 0,  and  u^(C)  < 0.  Thus, 

u < 0,  on  BC  . 


(3.9) 


he  functions  and  -uy  are  harmonic  conjuuates  in  fi  and  u is  bounded. 


V 


Applying  Lemma  2.1  and  noting  (3.2)  as  well  as  (3.7)  through  (3.10),  we  conclude  that 

Ux  = Ux  0 in  A • (3.10) 

Thus,  u^  ■ 0 in  il  since  otherwise  the  strong 


maximum  principle  (Courant  and  Hilbert 
(b>62,  p.  3201)  implies  that  Uj{  C in  fl.  which  is  not  possible.  Hence 

|grad  u|'  * u‘  ' 0 in  i)  . (3.11) 

Remark  1^1 

l-  decreases  monotonically  on  ABC  and  increases  monotonical ly  on  CPE.  The  tact 

that  |grad  u|  * 0 in  C is  also  a consequence  of  a result  of  Walsh  11^50,  p.  318, 

last  paragraph)  on  the  critical  points  of  harmonic  functions.  □ 

Since  v o in  the  function 

tTn(-U  4 iu  ) = --  Cnfu*'  *■  * i arc  tan(-U  V ) , ( 3 * 1 J > 

x y x y y x 

is  regular  in  il  and  has  bounded  imaginary  part  in  From  lemma  2.1  and  (3.2) 

through  (3.4)  we  conclude  that 

2 2 2 

Q - U ♦ li  ‘ % I '4 , in  i?  . (Ill) 

x y v . i 


Finally,  since  in  £>  is  harmonic  in  and  Q attains  its  minimum  on  F 


the  Hopf 


Principle  implies  that 

i.  P = I |2  < o,  on  r . 

We  now  use  Lemma  2.2  to  express  3Q/3 n in  terms  of  the  boundary  data  for  u 

on  T and  the  curvature  < of  T.  Direct  computation  shows  that,  on  T, 

U = n/2;  U = u = n/2  ; 
t s 

U = u - 7 n ^ ; n f = t E/2  , 

n n 2 1 3x  2 2 3y 

Utt  = Uss  + KUn  = '^/2  + K^/2  = ^ • 


Unt  = °ns  - KUs  = 5/2  - KT1/2  = -Ktl  ' 


0 = -U.  . 

nn  tt 


-kC  • 


2 2 2 

Thus,  since  Q = U + U , 
n t 


QQ  = U U + U U 
n n nn  t nt 


= -<c2/2  - <n2/2 


= -K/2  . 

Together,  (3.14)  and  (3.15)  show  that  < > 0 so  that  f is  convex. 


(3.14) 


(3.15) 
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Append i x : Remarks  on  the  doner all  zed  Max imum  I'tinciplo 

It  is  natural  to  try  to  weaken  tire  assumptions  of  Lemma  2.1  be,  .ruse  this  would 
sim;  lify  applications  to  free  boundary  problems.  For  example,  in  this  papet  flu 
assumption  that  ll  is  bounded  (see  (1.14)1  was  only  needed  so  that  lemma  2.1  could 

y 

be  applied.  The  proof  that  ll^  is  bounded  (Friedman  114761)  is  not  straightforward, 
ami  it  would  be  nice  if  it  could  be  avoided. 

The  requirement  in  lemma  2.1  that  v be  bounded  ean  be  replaced  by  the  weaker 
condition  that  v(u>(5))  r h^,  but  this  is  not  a very  convenient  hypothesis  to  check. 
Alternatively,  one  could  require  that  u + iv  r Ej  (ft)  (Goluzin  11469,  p.  4381, 
Priwalow  11956,  p.  1881). 

In  applications  we  will  usually  know  that 

2.2 


//  ( u + v2)dxdy  ' •»  , 

ft 

so  that  u and  v both  belong  to  the  space  L^Hfft)  which  consists  (Hille  11462, 


(A.l) 


p.  325))  of  functions 


which  are  harmonic  in  ft  and  satisfy 

//  ^“dxdy  < • . 

ft 


(A. 2) 


tt  is  tempting  to  conjecture  that  if  u,v  c L^H(ft)  then  u,v  r h,,  but  this 
is  not  true  as  the  following  example  shows.  Let  ft  * P,  so  that  c - z.  Set 


u(z)  ♦ iv(z)  • J z"/.tT  . (A.l) 

n-1 

Expanding  u and  v in  Fourier  series  and  using  the  orthogonality  relations  ot 

sin  nS  and  cos  n(l  we  see  that 

2*  . 2n  ", 

/ u2(r,e)dP  = / v (r,0)d0  » £ t ' n/n  . IA.4) 

0 0 n-1 


Thus  u,v  r but  u,v  f h^.  Furthermore,  usinci  a result  given  by  Tit « lunar sh 

119  39,  n.  lf>l,  Problem  1^1,  we  see  that 


Gohnncf  (1957)  (see  also  Tsuji  Il9r>9,  p.  186))  usrs  the  rn.ixim.il  theorem  of  H.iiily 
and  Littlewood  to  prove  that  it  u c L;>H(D)  then,  for  almost  all  9, 
u { z ) *o(li  - |z|)  ^ as  z * in  any  fixtni  Stole  domain  with  vertex  < . 

Unfortunately,  this  result  is  not  quite  strong  enough  to  prove  that  u is  bounded, 
even  if  we  assume  that  u is  bounded  on  except  at  z - 1.  otherwise,  by  mapping 

D onto  the  right  half  plane  we  could  conclude  that  a harmonic  function  on  the*  tight 
half  plane  which  has  bounded  limits  on  the  imaginary  axis,  and  bounded  growth  ft  ear 
the  teal  axis,  must  be  bounded.  This  would  constitute  a substantial  st lengthen i ng 
of  the  Phragmen-Lindelof  Principle  (Plotter  and  Weinberger  11967,  p.  94)). 

In  conclusion,  we  observe  that  the  behaviour  of  solutions  of  elliptic  equations 
near  corners  has  been  considered  by  Oddson  11968),  Kondrat  ' ev  1 196  7),  Miller  1 1'*07 , 
19711  and  (Ji isvard  (1969).  These  results  are  not  immediately  applicable  to  the 
present  problem  because  they  require  knowledge  of  the  behaviour  of  the  free  boundary  V 
near  its  endpoints.  Of  course,  once  it  is  known  that  T is  differentiable  at  its 
endpoints  then  this  work  will  yield  useful  information  about  the  behaviour  of  the 
solution  can  be  obtained. 
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